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Abstract 

We study the problem of selecting between different mechanisms of supersymmetry 
breaking in the Minimal Supersymmetric Standard Model using current data. We 
evaluate the Bayesian evidence of four supersymmetry breaking scenarios: mSUGRA, 
mGMSB, mAMSB and moduli mediation. The results show a strong dependence on 
the dark matter assumption. Using the inferred cosmological relic density as an upper 
bound, minimal anomaly mediation is at least moderately favoured over the CMSSM. 
Our fits also indicate that evidence for a positive sign of the ^ parameter is moderate at 
best. We present constraints on the anomaly and gauge mediated parameter spaces and 
some previously unexplored aspects of the dark matter phenomenology of the moduli 
mediation scenario. We use sparticle searches, indirect observables and dark matter 
observables in the global fit and quantify robustness with respect to prior choice. We 
quantify how much information is contained within each constraint. 
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1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM)[T], [2] with R-parity can solve the 
hierarchy problem and provide a candidate dark matter (DM) particle. The absence of 
observed supersymmetric (SUSY) particles in collider experiments to date implies that if 
nature is ever supersymmetric then it must be broken. The source of SUSY breaking is a 
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priori unknown in a bottom-up approach to particle physics phenomenology. As spontaneous 
SUSY breaking in the MSSM is not viable [2], one is led to consider SUSY breaking in a 
hidden sector of the theory communicated via a messenger sector to the MSSM. The choice 
of messenger sector leaves an imprint on the pattern of SUSY breaking, and therefore on 
the expected SUSY phenomenology. There are currently many viable soft-SUSY breaking 
schemes with some common mediation mechanisms including gauge, gravity, anomaly and 
moduli mediation. There are examples of each class of model which only have a few param- 
eters and can be somewhat constrained by current cosmological and indirect collider data. 
Sparticles appear in loop contributions to electroweak and 5-physics observables, affecting 
their values indirectly. In this paper, we ask the question: is there sufficient power in such 
combined data to favour one simple model over the other? Such model selection is likely to 
become even more interesting and important if signals compatible with SUSY are found in 
the Large Hadron Collider (LHC) experiments. 

To date, global Bayesian fits have been performed to the constrained MSSM (CMSSM) [3l 
SElElElElEllinillllllS], large volume string compactification models (LVS) ^LSj, the 
phenomenological MSSM (pMSSM) [T3] and the Non- Universal Higgs model [13]. In any 
Bayesian fit it is essential to check for robustness by choosing widely different but reasonable 
prior distributions of the model parameters. For models with more parameters the prior 
dependence becomes greater, but more precise and direct data will reduce such dependence. 
The pMSSM has twenty extra-Standard Model (SM) parameters and displays a large prior- 
dependence in the fit result^. Even the CMSSM, with only four extra-SM parameters, shows 
significant prior dependence and so many aspects of the fits are not robust. The parameter- 
space fits of the CMSSM and the pMSSM could only be robust with strong enough direct 
and precise data. The LVS fit, which has two extra-SM parameters, shows approximate 
prior independence [13]. Here, we shall perform fits to two additional models: minimal 
anomaly mediated SUSY breaking (mAMSB) and minimal gauge mediated SUSY breaking 
(mCMSB), each of which have the intermediate number of three parameters. Thus, we ask if 
the current available data are powerful enough to robustly constrain these models with three 
parameters additional to the SM as well as ask if it is powerful enough to disfavour any of the 
models over the others. The MSSM also contains a fi parameter in the Higgs superpotential, 
whose magnitude is constrained by the Z boson mass but whose sign is unknown. In the 
past, MSSM Bayesian model selection has only focused on the extent to which /i > is 
favoured over /i < [6l [131 [El [IB]. Thus, our work extends the use of Bayesian model 
selection to cover different SUSY breaking scenarios. 

Aside from the Bayesian fits, there have also been some global profile likelihood (equiv- 
alent to minimizing x^) analyses of the CMSSM [T] [TUl [13 |TB|. Prior dependence does not 
appear in frequentist interpretations of fits and so one cannot be sure whether they are ro- 
bust. The Bayesian analyses indicate that such fits are not yet robust with current data. 
A recent analysis [19] compared the CMSSM, mCMSB and mAMSB using electroweak 
and i?-physics observables while omitting the dark matter constraint. The analysis found 
that mAMSB was slightly preferred over the two other models by Ax^ = 2 and that light 

^When we discuss the number of parameters additional to the SM in a model, we refer explicitly to 
continuous parameters only. 
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SUSY was slightly favoured by each model. We extend this work in several ways, the most 
important being that we include dark matter as a constraint. We shall show that the dark 
matter constraint contains more information content than the other observables. By per- 
forming a Bayesian fit with several priors we are able to check robustness of the fit, and, 
unlike Ref. [19], we incorporate the uncertainties of important SM parameters into our fit. 

The parameters of the CMSSM are a flavour blind SUSY breaking scalar mass mo, a 
common gaugino mass M1/2, a flavour blind SUSY breaking scalar trilinear coupling Aq and 
tan P, the ratio of the MSSM Higgs vacuum expectation values (VEVs). Below a grand uni- 
fication theory (GUT) scale of Mqut ~ 2 x 10^^ GeV, the SUSY breaking terms of different 
flavours evolve separately to the weak scale. In anomaly mediated SUSY breaking (AMSB) 
[20J SUSY-breaking is communicated to the visible sector via the super- Weyl anomaly. In its 
original manifestation, pure anomaly mediation suffers from negative slepton mass squared 
parameters, signalling a scalar potential minimum inconsistent with a massless photon. Min- 
imal AMSB (mAMSB) assumes the existence of an additional contribution to scalar masses 
mo at Mgut giving it a total of three parameters: the VEV of the auxiliary field in the 
supergravity multiplet representing the overall sparticles mass scale, ruaux^ tt^o and tan/3. 
As advertised above, minimal gauge mediated SUSY breaking (mGMSB) [21] also has three 
continuous parameters: the overall messenger mass scale, Mmessi a visible sector soft SUSY- 
breaking mass scale, A and tan (3. It also contains an additional discrete parameter, namely 
Nmessi the number of SU(5) 5 © 5 representations of mediating fields. The example of a 
moduli mediated model which we consider is the Large Volume Scenario (LVS) derived in 
the context of II B flux compactification [22l [231 [2H [25l [13] , whose two extra-SM parameters 
are an overall SUSY breaking mass scale M3/2 and tan/5. 

The CMSSM is the phenomenologically most studied model in the literature, and so it 
is useful as a benchmark as to how well other models fare in the fits. It is not clear without 
further model building how the GUT-scale flavour blind structure of the CMSSM emerges, 
however. The main motivation for mAMSB and mGMSB is that the flavour problem is 
solved, since gauge interactions in mGMSB necessarily lead to flavour blind soft terms pi] 
at the messenger scale, and the mAMSB soft-terms are dominated by flavour diagonal pieces 
proportional to the gauge couplings [26] . The main motivation for LVS is that it results from 
a string compactification scenario which has moduli stabilization. 

The rest of the paper proceeds as follows: in Section [2] we show the technical aspects 
of our analysis, including a discussion of model parameters, prior probability distribution 
functions (PDFs), observables and the calculation of the likelihood. Section [3] presents an 
analysis of SUSY breaking model selection including the Bayesian evidence results on the 
preference for the sign of /x, the model selection and the dependence of our results on the 
dark matter constraint. Sections H] and [S] detail the phenomenologically viable regions of the 
parameter spaces of mGMSB and mAMSB, while Section [S] discusses the effects of individual 
observables and presents the best-fit points obtained, concluding with a discussion of the 
interplay between the anomalous magnetic moment of the muon and the rare branching 
ratio BR{B ^s7) in the different models. Finally, we discuss the implications of our fits 
for the direct detection of dark matter in Appendix \^ 
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1.1 Bayesian inference 

Problems in data analysis generally divide into two categories: parameter estimation and 
model selection. In parameter estimation problems one is interested in making inferences 
about the parameters of a given model using the available data and any other prior informa- 
tion. Model selection problems are concerned with distinguishing between different theories 
describing a given phenomenon. For instance, in the case of the CMSSM, one would like to 
know whether there is sufficient evidence in the data to rule out the fi < branch. Bayesian 
inference provides a consistent approach to model selection as well as to the estimation of a 
set of parameters m in a model (or hypothesis) H for the data d. It can also be shown that 
Bayesian inference is the unique consistent generalisation of the Boolean algebra [27] . 

Assuming some model hypothesis H, Bayesian statistics helps update some PDF |j(m|i7) 
of model parameters m with data. The prior encodes our knowledge or prejudices about the 
parameters. Since p(m|if) is a PDF in m, J p{m\H)dm = 1, which defines a normalization 
of the prior. One talks of priors being 'flat' in some parameters, but care must be taken 
to refer to the measure of such parameters. A prior that is fiat between some ranges in a 
parameter mi will not be fiat in a parameter x = log mi, for example. The impact of the 
data is encoded in the likelihood, or the PDF of obtaining data set d from model point nr. 
p{d\m,H) = C{m). The likelihood is a function of x^? i-e. a statistical measure of how well 
the data are fit by the model point. The desired quantity is the PDF of the model parameters 
m given some observed data d assuming hypothesis H: p{rn\d,H). Bayes' theorem states 
that 

/ I , o-N P{d\m, H)p{m\H) 
pim\d,H) = , (1) 

where p{d\H) = Z is the Bayesian evidence, the probability density of observing data set d 
integrated over all model parameter space. The Bayesian evidence is given by: 



where the integral is over dimensions of the parameter space m. Since the Bayesian 
evidence is independent of the model parameter values m, it is usually ignored in parameter 
estimation problems and posterior inferences are obtained by exploring the unnormalized 
posterior using standard Markov Chain Monte Carlo sampling methods. 

In order to select between two models Hq and Hi one needs to compare their respective 
posterior probabilities given the observed data set d, as follows: 



where p{Hi)/p{Hq) is the prior probability ratio for the two models, which can often be set 
to unity but occasionally requires further consideration (see e.g. [28t [29] for cases where the 
prior probability ratios should not be set to unity). It can be seen from Eq. [3] that Bayesian 
model selection revolves around the evaluation of the Bayesian evidence. As the average of 
likelihood over the prior, the evidence automatically implements Occam's razor. A theory 




(2) 



p{Hi\d) ^ p{d\Hi)p{Hi) ^ Zi pjHi) 
p{HM) pid\Ho)pm ZopiHoY 



(3) 
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|AlogZ| 


Odds 


Probability 


Remark 


< 1.0 


<3:1 


< 0.750 


Inconclusive 


1.0 


~ 3 : 1 


0.750 


Weak Evidence 


2.5 


~ 12 : 1 


0.923 


Moderate Evidence 


5.0 


~ 150 : 1 


0.993 


Strong Evidence 



Table 1: The scale we use for the interpretation of model probabilities. Here the 'log' 
represents the natural logarithm. 

with less parameters has a higher prior density since it integrates to 1 over the whole space. 
Thus, there is an a priori preference for less parameters, unless the data strongly require 
there be more. 

Unfortunately, evaluation of Bayesian evidence involves the multi-dimensional integral 
in Eq. [2] and thus presents a challenging numerical task. Standard techniques like ther- 
modynamic integration [30] are extremely computationally expensive which makes evidence 
evaluation typically at least an order of magnitude more costly than parameter estimation. 
Some fast approximate methods have been used for evidence evaluation, such as treating the 
posterior as a multi-variate Gaussian centered at its peak (see e.g. but this approxima- 

tion is clearly a poor one for highly non-Gaussian and multi-modal posteriors, such as those 
we sample here. Bridge sampling [321 [33l [M] allows the evaluation of the ratio of Bayesian 
evidence of two models, but can yield rather imprecise results. The problem can however 
be tackled by bank sampling but it is not yet clear how precisely bank sampling can 
calculate the evidence ratio. Various alternative information criteria for model selection are 
discussed by [36], but the evidence remains our preferred method. 

The nested sampling approach, introduced in [37], is a Monte Carlo method targeted 
at the efficient calculation of the evidence, but also produces posterior inferences as a by- 
product. [381 ES] built on this nested sampling framework and introduced the MultiNest 
algorithm which is efficient in sampling from multi-modal posteriors that exhibit curving 
degeneracies. MultiNest produces posterior samples and calculates the evidence and its 
uncertainty. This technique has greatly reduced the computational cost of model selection 
and the exploration of highly degenerate multi-modal posterior distributions. We employ 
this technique in this paper. 

The natural logarithm of the ratio of posterior model probabilities provides a useful guide 
to what constitutes a significant difference between two models: 



A log Z = log 



p{Ho\d) 



log 



(4) 



We summarize the convention we use in this paper in Table [H 
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2 The Analysis 



2.1 Choice of Prior Probability Distributions 

While for parameter estimation, the priors become irrelevant once the data are powerful 
enough, for model selection the dependence on priors always remains (although with more 
informative data the degree of dependence on the priors is expected to decrease, see e.g. |1Q]); 
indeed this explicit dependence on priors is one of the most attractive features of Bayesian 
model selection. Priors should ideally represent one's state of knowledge before obtaining 
the data. Rather than seeking a unique 'right' prior, one should check the robustness of 
conclusions under reasonable variation of the priors. Such a sensitivity analysis is required 
to ensure that the resulting model comparison is not overly dependent on a particular choice 
of prior and the associated metric in parameter space, which controls the value of the integral 
involved in the computation of the Bayesian evidence (for some relevant cautionary notes on 
the subject see ^T]). 

We have considered three different prior probability density functions (PDFs) in this 
analysis. The first is the standard "linear prior" where 

p(mi) = p{m2) 

for mi^2 two different points in the parameter space of one of the models under consideration. 
In particular, the linear priors are flat in the ratio of the two MSSM Higgs vacuum expectation 
values (VEVs), tan/5. It is important to realize that a prior which is fiat in one set of 
parameters m is not necessarily fiat in a different set of parameters m', say. The two priors 
will be related by a Jacobian factor such that 



p{m) 



dm' 



dm 



p{m'). (5) 



One may consider a more fundamental set of parameters to be those involving the quantities 
that actually appear in the Lagrangian before spontaneous symmetry breaking. Such a set 
would be fiat in the Higgs parameters B and /x, but not in tan p. "Natural priors" [7j are 
priors which are fiat in B and /i. The relationship between the two sets of parameters is 
given by the conditions of electroweak symmetry breaking 

„ sin2/? . , 5 9. , , mjr^ - mjjtan^ /3 Ml , , 

/iE = [m\^ + m\^ + 2/i2) and /i^ = \^^2 f_^ - -f^ (6) 

which lead to the Jacobian factor [7] 

B tan^ p - 1 



/itan /? tan^ /5 + 1 



(7) 



One may go further to examine the dependence upon the Higgs VEVs, Yukawa couplings 
and renormalisation group effects in more sophisticated priors [12]. We instead focus on a 
few reasonable but sufficiently different priors in order to check the prior independence of 
any inference we make. 
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CMSSM 


mAMSB 


50 GeV < mo < 4 TeV 
50 GeV < mi/2 < 2 TeV 
-4 TeV < Ao < 4 TeV 


50 GeV < mo < 4 TeV 
20 TeV < m3/2 < 200 TeV 


mGMSB 


LVS 


10^ GeV < A < 10^ GeV 
1.01 < M^e.s/A < 2 X 10^ 


50 GeV < mo < 2 TeV 



Table 2: Ranges for the parameters in mGMSB and the Large Volume Scenario. In mGMSB 
we also vary the discrete parameter Nmess between 1 and 8. For all models, 2 < tan/? < 62. 

It will also be the case in this analysis that we wish to estimate the scale of a parameter, 
rather than its exact value. This is the case for the messenger scale Mmess in gauge mediation, 
for example. In this case the appropriate prior to use is fiat in log(mi), and the relevant 
Jacobian factor is proportional to m|f^. We call this the "Jeffreys prior" or simply the "log 
prior" . 

2.2 Parameters and ranges 

Before proceeding we specify the parameter ranges over which we sample for the different 
models. Firstly, we consider both signs of fi in our analysis for all models. The ranges over 
which we vary the continuous model parameters are shown in Table El tan j3 is bounded from 
below by 2, values lower than this are in contravention of LEP2 Higgs searches, and from 
above by 62, since such large values lead to non-perturbative Yukawa couplings below the 
GUT scale and calculability is lost. In mGMSB the discrete parameter Nmess, the number of 
messenger multiplets, is varied between 1 and 8. Higher values of Nmess lead to problems with 
perturbativity of gauge interactions at the GUT scale |21j. In the CMSSM the unification 
scale is the standard GUT scale mcur ~ 2 x lO^^GeV, while for the LVS the soft terms are 
defined at the intermediate string scale m^ ~ lO^^GeV as in |13] . 

2.3 The likelihood 

The calculation of the likelihood follows [7] with updated data and additional observables. 
The constraints we use are all shown in Table [3] along with their respective experimental 
and theoretical sources. We treat the measurements Di of the observables as independent 
and we have been careful to eliminate possible correlations between them. We also assume 
Gaussian errors on all measurements except where explicitly mentioned. The log likelihood 
of a prediction pi of an observable i is given by 

logA = -^-^log(27r)-log(a,:) (8) 
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[67, 68J 


[69, 63J 


Rbtu 


1.28 ±0.40 


[61, 62, 67J 


[63] 


Ao- 


U.UOi_Q Q25 


[61, 62J 


[551 [701 [71] 


Rm 


1.004 ±0.007 


^l62j 


m 




172.4 ± 1.2 




m 




4.20 ±0.07 








91.1876 ±0.0023 




[551 [591 EO] 




0.1172 ±0.002 




m 




127.918 ±0.018 







Table 3: Experimental constraints, showing the observables, the constraints applied and the 
source of the theoretical and experimental values and errors. The first section shows general 
observables, the second electroweak observables, the third section B-physics constraints and 
the final section contains SM parameter constraints. The 95% confidence level (CL) direct 
sparticle search constraints from ^55j were also applied where relevant. 
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0.0343 0.1143 0.1943 

Figure 1: Depiction of our likelihood constraint on the predicted value of Qdm^^ due to the 
lightest neutralino, compared with a simple Gaussian with WMAP5 central value and a la 
uncertainty of 0.02. 



where Xi = l£L^iL^ q. jg the central experimental value and ai is the standard deviation 
incorporating both experimental and theoretical uncertainties. All likelihoods except for 
those of the DM relic density, the Higgs mass and BR{Bs — ^ /U^/U^) are calculated in this 
way. Details of the constraints applied on the last two exceptions mentioned above can 
be found in [TB]. The £j considered here are independent, and are multiplied together to 
calculate the overall likelihood. 

For f^DM^^, we use two different forms of constraint. Assuming that the DM relic density 
is composed entirely of the lightest SUSY particle (LSP), one obtains a Gaussian £dm with 
the WMAP5 central value. We call this the 'symmetric £dm'- It is of course possible that 
the neutralino does not make up the entirety of the DM in the universe. For instance, axions 
or stable moduli could make up the additional component. Indeed, in the Large Volume 
Scenario, the properties of the lightest Kahler modulus as a DM constituent have been 
discussed in [75]. In mAMSB, the wino co-annihilation is so efficient that the relic density is 
generically far below the WMAP value, so that such an extra component is cosmologically 
necessary. This motivates the use of the 'asymmetric £dm' which is given as: 

if X < c 
if X > c. 



exp 



2s2 



where c and s are the mean and 1-a error values of the Gaussian £dm- We take c = 0.1143 
that results from a fit to WMAP5, baryon acoustic oscillations and Type la supernovae 
observations. We also take s = 0.02 in order to incorporate an estimate of higher order 
uncertainties in its prediction [TS] . A diagram of the resulting likelihood penalty is displayed 
in Fig. [H 

In mGMSB the gravitino is the LSP. The gravitino mass 771^/2 depends on the scale 
of SUSY breaking and is usually in the range of a few eV up to a few hundred keV. This 
property has often been cited as an attractive feature of mGMSB theories, as the low mass of 
the gravitino indicates that gravity mediated effects are much smaller than gauge mediated 
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effects, so that ffavour changing neutral currents are naturally suppressed. The gravitino 
contribution to the relic density is given by |77| 




100 



(10) 



keV [g^iTf) 



where g*{Tf) is the number of massless degrees of freedom at the gravitino freeze-out temper- 
ature, and for SUSY models is in the range 100-200 [21]. We see from Eq. [10] that achieving 
the WMAP value of the relic density requires a gravitino at the lower end of the mass range, 
in particular less than 1 keV. Constraints on structure formation and WMAP [78] are now 
strong enough, however, to rule out the resultant warm DM. Heavier gravitinos which evade 
this bound would be over abundant compared to observation. Although it may be possible to 
dilute the gravitino relic density by late entropy production, we wish to keep our analysis as 
general as possible without confining ourselves to a specific change in the physics of the early 
universe. Accordingly, in our analysis we do not impose any DM constraint on mGMSB. 

The combined log likelihood is the sum of the individual likelihoods. 



To calculate the MSSM spectrum we use Sof tsusy2 .0.18 [79j which calculates the spectrum 
of the CMSSM, mAMSB and mGMSB. By modifying the unification scale from itigut to 
f^string ~ lO^^GcV and the gauge coupling boundary conditions Softsusy can also provide 
the spectrum in the LVS case. Parameter space points which violate the current sparticle 
exclusion bounds of [SS], do not break electroweak symmetry correctly or have tachyonic 
sparticles are assigned zero likelihood. Points which have a charged LSP are rejecteqj. If 
a point survives the cuts above, it is passed via the SUSY Les Houches Accord |81] to 
microMEGAS2.2 jUl UHl HI, |50] , SuperIso2.3 |6ll[62] and SusyPOPE [S]. From microMEGAS 
we obtain the DM relic density, the rare branching ratio BR{Bs — > fi~^fi~), the SUSY compo- 
nent 5a^ of the anomalous magnetic moment of the muon {g — 2)^ and DM direct detection 
rates. To the one-loop value of calculated by microMEGAS, we add the logarithmic piece 
of the quantum electro-dynamics 2-loop calculation, the 2-loop stop-Higgs and chargino- 
stop/bottom pieces and the recently discovered two-loop effect due to a shift in the muon 
Yukawa coupling and proportional to tan^ (3 [521 [531 [SI] • From SuperIso2 . 3 the branch- 
ing ratios BR{B Xgj), BR{B — >• Dtp), the quantities RaMs^ Rm, Rbtu and the isospin 
asymmetry Aq- are obtainecif|. SusyPOPE provides all of the observables listed in the second 
part of Table [31 

We present our data in the form of marginal posteriors in one and two dimensions, 
where the unseen dimensions have been integrated over. To do this integration numerically 

^Due to the small ncutralino-chargino splitting in niAMSB we have been careful to reject any points 
that would violate the long-lived charged stable particle bounds from Tevatron, which requires Am = 
— m^o > 50 MeV. In fact, we find that this bound does not constrain the mAMSB parameter space 
since mAMSB predicts larger splittings [50] , 

■^We note that in the process of preparing this paper and after our fits were performed a new version 
of SusyBSG[82j appeared. This more accurate calculation could result in a change in our BR{B Xg^) 
prediction of up to 0.13 x 10^"*. 




(11) 
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we divide the range of the parameter in which we are interested into a series of bins, and 
then count the number of samples which fall into each bin. We use 60 bins per parameter, 
which is a trade-off between parameter resolution and unwanted statistical noise from finite 
sampling effects. We shall also discuss the profile likelihood at various stages. The profile 
likelihood is often shown in one dimension, and plots the maximum of the likelihoods of the 
samples that fall in each bin. This statistic is equivalent to plotting minimum in various 
directions of parameter space, and shows where the data are fitted best, but does not take 
into account volume effects of the less best-fit points. Everywhere in this paper, we take 
best-fit to mean the highest likelihood (or equivalently, lowest x^)- MultiNest is optimized 
for calculating Bayesian evidence values, and also typically gives a reasonable sampling for 
posterior evaluation. However, the profile likelihood is extracted by MultiNest with a lot 
of noise. It was found by Ref. [10] to depend upon the prior, which it should not do for 
a large enough sampling. Therefore, here we do not explicitly include plots of the profile 
likelihood, preferring instead to mention its behaviour at various points in the text. 

We have taken care that all the two-sided constraints applied in the likelihood are nor- 
malised to 1. For one-sided constraints such as sparticle exclusion limits, the likelihood is 
zero below the 95% confidence limit, and 1 above that limit. The likelihoods we apply for 
BR{Bs Ai+yU^) and m/j are more complicated, but they both asymptote to 1 in the limit 
of zero branching ratio and large Higgs mass respectively. While we focus on the statistically 
meaningful difference between evidences in our fits, by specifying the normalisation of the 
likelihood and providing also the absolute values of the evidence in our tables of results we 
hope to encourage other groups, by following our calculation, to directly compare Bayesian 
evidences calculated for other models with the results presented herein. 

3 Model Selection 

In this section we discuss whether current data prefer one model over the others in the set 
mAMSB, mGMSB, the CMSSM and LVS. We also present our results on quantification of the 
preference of the fits for yU > 0. In both these cases, we do the analysis with 'symmetric £dm' 
as well as 'asymmetric £dm' (see Section l273l) . Strictly speaking, a calculation of the Bayesian 
evidence with asymmetric dark matter constraints is inherently unfair unless the relic density 
of the additional DM component is included, with an associated reduction in the evidence 
due to the Occam's razor effect. We shall here make model selections between models 
where the asymmetric likelihood penalty has been applied to all of the models in question. 
Thus, although the Occam's razor effect is neglected, it should affect all of the models being 
compared in approximately the same way and so the comparison should remain unaffected. 
We also assume that the prior probabilities of each model are equal and independent of the 
sign of n- 
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symmetric £dm 


asymmetric £dm 


Model/Prior 


linear 


log 


natural 


linear 


log 


natural 


CMSSM 
mAMSB 
LVS 


8.0 ±0.1 
0.4 ±0.1 
8.7±0.1 


7.9 ±0.1 
0.6 ±0.1 
8.9 ±0.1 


10.3 ±0.1 
0.0 ±0.1 
11.8 ±0.1 


0.0 ±0.1 
5.1 ±0.1 
2.9 ±0.1 


1.0 ±0.1 
6.0 ±0.1 
3.0 ±0.1 


1.3 ±0.1 

5.0 ±0.1 

3.1 ±0.1 



Table 4: log evidences (Alog^) for mAMSB, LVS and the CMSSM for both signs of ^. 
Symmetric Cbm labels the assumption that the DM relic density is composed entirely of 
the LSP and asymmetric Cdm labels the assumption that the LSP forms only a part of the 
DM relic density. The log evidence of the natural prior mAMSB, logZg = 67.3 and the 
log evidence of the linear prior CMSSM, 76.7 have been subtracted from all entries in the 
symmetric £dm and asymmetric £dm respectively. 



3.1 A preference for mAMSB, mGMSB, the CMSSM or LVS? 

We now discuss to what extent the mAMSB, mGMSB, the CMSSM, and LVS model are 
preferred over one another by the current data. In this subsection, we marginalise over the 
sign of /i. 

The log evidence values (log Z) for mAMSB, LVS and the CMSSM are the most important 
results of this paper and they are shown in Table |H mAMSB is strongly favoured over the 
CMSSM (A log Z > 5) and moderately favoured over LVS (log AZ > 2) for asymmetric Cbm, 
a result which is approximately prior independent. However, mAMSB is almost decisively 
ruled out for symmetric Cbm (A logZ < —7). Although mAMSB with a purely thermal relic 
density is decisively disfavoured we have not taken into account the non-thermal component 
of the relic density due to decays of the heavy gravitino in mAMSB which for, some values 
of the inflationary reheating temperature, could saturate the WMAP bounds. As mentioned 
above, in mAMSB, the degeneracy between the LSP and next-to-LSP (NLSP) chargino 
causes extremely efficient co-annihilation of the LSP to near zero relic density over the whole 
parameter space. Thus, the entire parameter space of mAMSB passes the asymmetric Cdm 
constraint and it is strongly preferred. On the other hand the LSP relic density in most of the 
mAMSB parameter space is much lower than the central value inferred by WMAP and hence 
mAMSB is almost ruled out if one assumes that the DM relic density is composed entirely 
of the LSP. There is a weak to moderate preference for LVS over the CMSSM depending on 
the prior and whether the symmetric or asymmetric form of Cbm is taken. 

To quantify the extent to which these results depend on the DM constraint, we calculate 
the Bayesian evidence ratios for mAMSB, mGMSB, LVS and the CMSSM without this 
constraint and list them in Table O It is clear from these evidence values that the results are 
completely inconclusive in the absence of the DM constraint and hence it can be concluded 
that it does indeed dominate our model selection results. Ref. [19] concluded on the basis 
of a minimisation with a subset of the observables included in the present paper and 
no DM constraint that mAMSB was slightly preferred over the other two models. Table [5] 
indicates that even such a weak inference cannot yet be made due to non-robustness of fits 
in mAMSB, since the preference for or against it depends upon the form of the prior. 
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Model/Prior 


linear 


log 


natural 


CMSSM 
niAMSB 
mGMSB 
LVS 


0.9 ±0.1 
1.1 ±0.1 
1.4 ±0.1 
1.1 ±0.1 


1.0 ±0.1 
1.2 ±0.1 
1.5 ±0.1 
0.8 ±0.1 


1.1 ±0.1 
2.9 ±0.1 
0.0 ±0.1 
1.5 ±0.1 



Table 5: log evidences (logZ) for mAMSB, mGMSB, LVS and the CMSSM for both signs 
of fi and without imposing the ^dm^^ constraint. A constant value (the log evidence of 
mGMSB using natural priors, 78.2) has been subtracted from all the log evidence values. 





symmetric £dm 


asymmetric £dm 


Model/Prior 


linear 


log 


natural 


linear 


log 


natural 


CMSSM 
mAMSB 
LVS 


1.2 ±0.1 
1.4 ±0.1 
3.2 ±0.1 


2.4 ±0.1 

2.5 ±0.1 
3.1 ±0.1 


0.4 ±0.1 
0.4 ±0.1 
2.6 ±0.1 


1.3 ±0.1 
1.9 ±0.1 
3.6 ±0.1 


2.4 ±0.1 
3.4 ±0.1 
3.9 ±0.1 


0.4 ±0.1 
0.6 ±0.1 
3.3 ±0.1 






No DM constraint 


mGMSB 




1.7 ±0.1 


2.2 ±0.1 


1.4 ±0.1 



Table 6: Difference in log evidences (A log 2) for > and /i < for the CMSSM, mAMSB 
and the LVS model. A positive value for a model indicates a preference for /i > 0. 'Symmetric 
£dm' labels the assumption that the DM relic density is composed entirely of the LSP and 
'asymmetric £dm' labels the assumption that the LSP forms only a part of the DM relic 
density. 

3.2 Selection of sign(;u) 

We summarize the amount of preference for /i > in mAMSB, the CMSSM and LVS in 
Table El We list the Bayes factors, A log Z in favour of /x > for these three models. Using 
the Jeffreys scale defined in Table [H it is evident that although there is a positive but not 
strong evidence in favour of /x > for all three models, the extent of the preference depends 
quite strongly on the priors used. Although there are small numerical differences due to the 
constraints being updated and different prior ranges taken, the CMSSM A log Z values lead 
to the same conclusions as previous determinations P [131 [HI [IS] : the preference for /i > 
is moderate at best, and quite prior dependent. This dependence on the prior is a clear sign 
that the data are not yet of sufficiently high quality to be able to distinguish between these 
models unambiguously and hence it is not justified to ignore the /x < branch in any of the 
models under consideration, despite the fact that ((7 — 2)^ favours positive fj, by around SAa. 
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A (GeV) ^^qS a (GeV) ^^^5 

Figure 2: 2D marginalised posterior PDFs for mGMSB in the log]^o(^mess)-A plane for (a) 
log priors and (b) natural priors with 68% and 95% Bayesian credibility intervals. Both plots 
are marginalized over both signs of fi. Green dots show the position of the best-fit point for 
each prior choice sample. 



Oh 




Figure 3: mGMSB ID marginalized posterior PDFs over (a) logiQ{Mmess) and (b) A for 
linear, log and natural priors. Both plots are marginalized over both signs of fi. Both here 
and in future one-dimensional marginalizations, the vertical axis has been normalized so that 
the maximum posterior PDF is one. 
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4 mGMSB Parameter Constraints 



We now discuss the constraints on mGMSB coming from the global fit. These have not 
appeared in the literature before, except for in Ref. [19], where a minimisation was 
performed, and so prior dependence was not quantified. Fig. [2] shows the marginalized joint 
posterior PDFs for mGMSB in the logiQ{Mmess)--^ plane for logarithmic and natural priors. 
We have also marginalized over the messenger index Nmess- The most striking feature of the 
plots shown is that while the 2D posterior with the log prior is flat in the logarithm of the 
messenger scale Mmess , for the natural priors there is a very strong peak at the highest values 
of Mmess- This is due to a volume effect. With a linear or natural prior, smaller values of 
Mmess occupy a smaller prior volume, which pushes the posterior for Mmess down. The log 
priors sample more from the region of low Mmess , due to the suppression from the log prior 
measure, and are in rough agreement with the profile likelihood. The excluded region at low 
values of the messenger scale in Fig. [2](a) is due to the requirement that Mmess > A, since if 
this is not the case the messengers are tachyonic. In order to understand why Mmess is only 
very weakly constrained by the data, consider the spectrum resulting from gauge mediation. 
To leading-log order, the gaugino masses at the scale Mmess are 

M-^XMmess) = hNmessA^^^^^^9{A/Mmess), (12) 

where hi = (5/3, 1, 1) and kiai (no summation) are all equal at the GUT scale, and a, are 
the gauge coupling constants. The messenger scale threshold function g can be found in [21] 
and g{x) 1 when A -C Mmess- This is the case throughout most of the mGMSB parameter 
space that we study. In this limit one can also approximate the scalar masses by 

m)iMmess) = 2Nmess g Qh^^^^A' - (13) 

where Cj are the quadratic Casimir operators of the relevant gauge groups [21] . The depen- 
dence on Mmess is somcwhat subtle. Increasing Mmess leads to different initial conditions 
for ai{Mmess) and also increases how far the renormalisation group equations (RGEs) must 
be evolved to reach the low scale. As pointed out in [83j these effects cancel each other 
out in the gaugino masses since the gaugino masses obey the same one-loop RGEs as the 
gauge couplings. There is however an effect on the scalar masses, leading to a decrease in 
the squark masses and an increase in the slepton masses at larger values of Mmess- This 
statement does not apply to the stop masses, which do depend somewhat on tan/? due to 
mixing. Therefore the mGMSB spectrum is fairly independent of Mmess- We remind the 
reader that for the mGMSB fits, no DM constraint is applied. 

Fig. [3]^b) shows the ID PDF posterior for the SUSY breaking scale A. All sets of priors 
have a strong peak at around 3 x lO'' GeV before falling away at higher values of A. Since 
A is directly related to the scale of the SUSY spectrum, high values of A lead to a heavier 
spectrum. Such a spectrum fits {g — 2)^ and BR{b 57) worse than a lighter spectrum, as 
we shall show in section 16. 3[ 

We next turn to tan/5, for which we present the posterior PDFs in Fig. ID^a). The 
posterior PDFs for the logarithmic prior have the same shape as the profile likelihood, with 
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tanp 



Figure 4: niGMSB ID marginalized posterior PDFs for (a) tan/5 and (b) the messenger 
index Nmess for log, linear and natural priors. 



both sharply rising from the allowed value of tan/3 = 2 to a most likely value just under 
20 and then gradually falling away as tan/3 increases, a feature common to both branches 
of /i. It is primarily the anomalous magnetic moment of the muon that is pulling the 
model towards high tan/?, where a significant SUSY component can fit the data. SUSY 
contributions to {g — 2)^ are dominated by one loop diagrams involving chargino-sneutrino 
or neutralino-smuon in the loop. These are approximately 

_ ^yvsY ^ + glF.M,)) (14) 

where is the muon mass and Fi^2 are positive definite functions of sparticle masses, 
which scale as l/Mg^gY the limit of heavy degenerate sparticles of mass Msusy- With 
natural priors the behaviour is similar, except that the peak occurs at a lower value of tan (3 
driven by the fact that natural priors prefer small tan/3, which can be seen by inspecting 
the Jacobian factor in Eq. [71 Fig. ID^b) shows the posterior PDFs for the discrete messenger 
index Nmess- The log priors prefer a large number of messenger multiplets, with the most 
likely value of Nmess clearly being 8. Since the log prior prefers light values of the neutralino 
and chargino masses, a good fit to the anomalous magnetic moment of the muon is also 
preferred by having relatively light scalars, so that the smuon and sneutrino diagrams may 
contribute significantly. From Eqs. [12] and [131 a slepton mass divided by a weak gaugino 
mass is proportional to Nmess and so high Nmess, i-e. light sleptons is preferred. This 
preference is less strong for the heavier spectra with the linear or natural priors, where it is 
more difficult to fit {g — 2)^. It also has to compete with volume effects coming from the 
LEP2 Higgs constraint, which prefers heavy stops and therefore low Nmess- 

To summarize, all of the mGMSB parameters show significant prior dependence, even 
when each is marginalized down to one dimension. We decline to present posterior PDFs 
of sparticle masses, since they show the same non-robustness with respect to changing the 
priors. 
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5 mAMSB Parameter Constraints 



In this section we discuss the effects of the observable constraints on the anomaly mediated 
SUSY breaking parameter space. All supergravity theories suffer from a Weyl (rescaling) 
anomaly which leads to soft breaking terms for the visible sector. These contributions are 
usually rendered negligible by the usual gravity mediated soft breaking terms. If the gravity 
mediated terms are suppressed, perhaps by a mechanism similar to that originally suggested 
in the brane model of [20] then the SUSY breaking is dominated by the anomaly mediated 
terms to leading order. In pure AMSB the slepton masses suffer from being tachyonic. One 
way this can be ameliorated while avoiding the SUSY flavour problem is by introducing a 
new universal mass parameter mo at the GUT scale which lifts the slepton masses. We 
define the GUT scale to be the scale of electroweak gauge unification. The GUT-scale soft 
breaking terms for the gauginos, scalars and trilinear couplings are then given by 

M, = ^m3/2, (15) 
9i 

^ ~\ ijt^'' ^ S^') ^ ^^^^ 

Ay = -^m3/2, (17) 

where j3g^ is the beta function of the z*^ SM gauge group, 7 is the anomalous dimension 
of the respective scalar wave-function and y and (3y are Yukawa couplings and Yukawa 
beta functions respectively. Eqs. [TMT7I constitute the "minimal AMSB" (mAMSB) SUSY 
breaking assumption. Its spectrum is specified by three continuous free parameters: the 
gravitino mass ^3/2, the parameter mo and tan/?, along with the sign of the /x parameter. 
Our mAMSB fits use the 'asymmetric' DM constraint shown in Fig. [H since the mAMSB 
neutralino relic density is too low to be the sole component of the dark matter. 

Fig. [5](a,c,e) shows the posterior PDFs for linear priors in the m^j2-TnQ, m3/2-tan/5 and 
mo-tan /9 planes while Fig. [ni^b,d,f) shows the same plots for natural priors. The dark matter 
constraint is asymmetric in both cases. There is clearly a strong prior dependence in all 
parameters. This is due to the pull of the natural priors to low values of tan (3 and to high 
values of -B//i, which occurs near the focus point at high mo |84i- The linear priors identify 
a different region of high probability at low mo and low m3/2 leading to a light sparticle 
spectrum capable of satisfying the constraints from {g — 2)^ and BR{B Xg'-f). Log prior 
posteriors are omitted for brevity, but they are similar to the linear prior results, except that 
the log priors have an even stronger preference for low mo and m3/2. 

There are three significant regions of parameter space which have been ruled out already. 
The first is the triangular region at low mo and moderate to high values of m3/2. This is 
ruled out by the direct search constraints on the sleptons. The large values of m3/2 and 
small mo lead to low (or even imaginary) values of the slepton masses in Eq. [TSl The 
disallowed strip at low values of m3/2 and all values of mo is due to gaugino masses being 
too low. The region defined by tan/5 > 50 is decisively disfavoured by a combination of 
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Figure 5: mAMSB marginalised 2D posterior PDFs in the mo-ms/2, mo-tan j3 and m3/2- 
tan jS planes for linear priors (left-hand side pictures) and natural priors (right-hand side 
pictures) with asymmetric Com- The green dots mark the best-fit points in each prior case 
sample. 68% and 95% Bayesian credibility regions are shown as the inner and outer contours 
respectively. 
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Figure 6: mAMSB marginalized ID posterior PDFs for (a) mo, (b) ^3/2 and (c) tan/? for 
different priors. Red (solid) lines represent linear priors, green (dashed) lines for log priors 
and blue (dotted) lines for natural priors. All plots are for asymmetric Com- 



related circumstances: tachyons, no radiative electroweak symmetry breaking and, as mo 
gets larger, a Higgs potential that is unbounded from below. 

Fig. [6] shows the ID posterior PDFs for (a) mo, (b) m3/2 and (c) tan/3 for all three priors 
used with red (solid line) representing the linear priors, green (dashed lines) the log priors 
and blue (dotted lines) the natural priors. The dark matter constraint is asymmetric in all 
cases. As expected the logarithmic priors show a distinct preference for low mo, with the 
linear priors expressing the same behaviour except for a longer and less suppressed tail to 
high mo- The natural priors have the opposite behaviour, indicating that the fits are strongly 
prior dominated. Despite the aforementioned statistical noise, profile likelihoods do indicate 
that the region of low mo favoured by the linear and log priors is a better fit to the data than 
the region favoured by the natural priors. Indeed, part of this peak is likely due to a volume 
effect due to the large probability mass at high mo- There is much better agreement between 
the three sets of priors for the posterior PDF for m3/2 with all priors exhibiting a peak in 
the likelihood at near 50 TeV falling off quite rapidly in the high m3/2 region. The posterior 
PDF for tan (3 shows two different peaks. One of these is due to the natural priors and is at 
low tan /3 ~ 8 due to the natural priors preference for low values of this variable. The other, 
broader peak is a better fit and occurs around tan/? ~ 22. This moderate value of tan/5 is 
preferred due to the statistical pull from observables such as [g — 2)^, BR{B —>■ Xs'j) and 
BR{Bs — > and values higher than this are disfavoured due to the SM constraints on 

these same observables as is discussed in more detail in [26] . 

To summarize, like mGMSB, all of the mAMSB parameters show significant prior de- 
pendence, even when each is marginalized down to one dimension. We decline to present 
posterior PDFs of sparticle masses, since they show the same non-robustness with respect 
to changing the priors. We turn now to a more detailed examination of the constraints from 
the observables. 
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6 Constraining Power and Statistical Pull of Observ- 
ables 



6.1 Information Content and Constraining Power 

In order to calculate the information gain in moving from a prior distribution to the posterior 
distribution, one can calculate the Kullback-Leibler (KL) [85] divergence (also called infor- 
mation divergence, information gain, or relative entropy) between the prior and the posterior 
distributions. Denoting the posterior distribution and prior distributions by p(m|o?, H) and 
p{rn\H) respectively, the KL divergence of observable i is defined as 

D,.-,^jpiMiHno,(^e^). (18) 

Dkl quantifies how much the prior PDF has been updated in its transformation to the 
posterior PDF. Using Eqs. fTSfTl the KL divergence of observable i is 



Pirn). 



dm, (19) 



where p{di\H) = J Ci{m)p{m\H)dm. Since the Standard Model observables are used in our 
fits as inputs, we consider them to be a part of the prior. The KL divergence of all of the data 
combined Dkl is easily obtained from the MultiNest algorithm which already calculates 
it to obtain the uncertainty on the evidence estimate (see [38] )• We defined the constraining 
power of observable i to be 

Cp, = (20) 

The KL divergence was also employed in [IQl to calculate the constraining power of observ- 
ables for the fi > branch of the CMSSM using linear and log priors and it was concluded 
that the information gain is dominated by the DM constraint which alone accounts for 
Cp = 0.8 — 0.95 depending on the prior. By examining the variation of the 95% CL regions 
in the CMSSM parameter space with and without the DM constraint, [IB] on the other hand 
concluded that the DM has little role in constraining the parameter space of the CMSSM. 
We update the analysis of [10] by including the additional electroweak observables returned 
by SusyPOPE and calculating Dkl for the observables for the CMSSM with both signs of fi. 
We list the Cp^ values in Table [7] for the CMSSM, mAMSB and LVS respectively. 

Table[7]shows that in the CMSSM, the DM constraint is dominant in constraining the pa- 
rameter space, as is familiar from current literature. Since the asymmetric Ldm rules out less 
of parameter space than the symmetric constraint, Cp „ is always smaller in the asym- 
metric case. However, i?— physics constraints dominated by BR{B — > X^j), electroweak 
constraints, the anomalous magnetic moment of the muon and the Higgs mass constraint 
all play significant roles too. By contrast, mAMSB, which predicts approximately zero DM 
relic density over its entire parameter space, is hardly constrained at all by DM constraints. 
In fact, none of the indirect constraints provide much constraining power except for the elec- 
troweak observables. Once one has taken the SM inputs into account, the indirect observables 
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symmetric £dm 


asymmetric £dm 


Observables/Prior 


linear 


log 


natural 


linear 


log 


natural 




CMSSM 




0.65 


0.57 


0.63 


0.56 


0.47 


0.56 


B-Physics 


0.44 


0.40 


0.41 


0.49 


0.45 


0.48 


BR{B ^ Xsi) 


0.43 


0.40 


0.40 


0.48 


0.44 


0.47 


Electroweak 


0.49 


0.36 


0.45 


0.53 


0.38 


0.55 




0.41 


0.32 


0.37 


0.44 


0.34 


0.48 


rrih 


0.42 


0.39 


0.39 


0.46 


0.42 


0.47 




mAMSB 




0.02 


0.01 


0.01 


0.01 


0.00 


0.00 


B-Physics 


0.03 


0.06 


0.03 


0.04 


0.07 


0.03 


BR{B ^ X,7) 


0.02 


0.05 


0.01 


0.03 


0.05 


0.02 


Electroweak 


0.12 


0.11 


0.14 


0.13 


0.11 


0.15 


Sa^ 


0.01 


0.05 


0.01 


0.02 


0.06 


0.02 


rrih 


0.02 


0.05 


0.04 


0.03 


0.05 


0.05 




LVS 




0.49 


0.45 


0.47 


0.26 


0.16 


0.23 


B-Physics 


0.19 


0.23 


0.18 


0.25 


0.27 


0.22 


BR{B ^ X,7) 


0.19 


0.22 


0.17 


0.24 


0.27 


0.23 


Electroweak 


0.14 


0.11 


0.13 


0.15 


0.11 


0.14 




0.09 


0.05 


0.08 


0.12 


0.06 


0.11 


rrih 


0.20 


0.20 


0.21 


0.25 


0.25 


0.26 



Table 7: Constraining power of different observables in the CMSSM, mAMSB and LVS 
marginalized over both signs of /i. 'B-Physics' includes the observables in the third section 
of Table [21 Electroweak observables include the ones in the second section of Table [3] plus 
and sin^ 9[jj:. Symmetric £dm denotes the assumption that the DM relic density is 
composed entirely of the LSP and asymmetric £dm codifies the assumption that the LSP 
forms only a part of the DM relic density. 
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add very little to the degree of constraint on parameter space. The LVS is most highly con- 
strained by DM, although less so than the CMSSM. This is because [13] a larger volume of 
the prior corresponds to a relic density compatible with the WMAP DM constraint compared 
to the CMSSM. All of the other indirect observables also help to constrain the LVS scenario 
in a non-trivial way. Defining Csm where SM = {as{Mz)^^ ,(y{Mz)'^^ ,mt,mf„mz}, we 
consistently find values Csm =0.8-0.9 for each of the three models investigated in Table [71 
regardless of the form of the DM constraint or prior. It is therefore essential to vary these 
input parameters and to use the available experimental data to constrain them. 

6.2 Statistical Pull of Observables 
6.2.1 mGMSB 

We now examine the pull exerted by the experimental constraints. Some of the more impor- 
tant of these are shown in Fig. [71 Considering first the posterior PDF for rrit, there is rough 
equality between the experimental constraint and that predicted by mGMSB, indicating that 
it is not being strongly pulled by other observables. 

The W boson mass is typically under predicted by about la as compared to current direct 
measurements of its mass. The fits show a weak correlation between rrit and mw such that 
a decrease in the experimental value of mw would also lead to a better fit for rrit. This effect 
is quite small however, as is the overall disagreement between theory and experiment. The 
PDFs for TTifo closely follows its experimental constraint and sin^ ^^^j is constrained to be close 
to its central value. These two observables are not shown, for brevity. In Fig. Wic) we show 
the logarithm of the branching ratio for the flavour changing decay Bg — > n^n^ . The 95% 
CL experimental upper limit from the Tevatron is marked with a black arrow at the right of 
the plot. If this branching ratio is measured to be non-zero in the near future, it will be far 
above the SM prediction and therefore provide a strong constraint on mGMSB, removing 
most of the currently available posterior density, which resides near the SM prediction of 

The 3.2(7 discrepancy between observation and the Standard Model theoretical value of 
the anomalous magnetic moment of the muon has received much attention from SUSY phe- 
nomenologists. In mGMSB the sign of the SUSY corrections to this quantity are correlated 
with the sign of /i. The ^ > branch thus fits the data much better than the fi < branch. 
For /i > the natural priors prefer small values of 6a^, since the SUSY contributions to 6a^ 
are proportional to tan (3, large values of which are suppressed by the natural prior. Log 
priors, while still exhibiting a preference for low values of 5a^, can provide a better fit. 

Fig. [71(e) shows that mGMSB predicts the most likely branching ratio of i? — > Xg'j to 
be between la less than the experimental central value. The stop-chargino contribution 
interferes (for positive values of M3 and At such as we have in mGMSB) opposite to the sign 
of fi with respect to the SM contribution. When > as preferred by 6a^, it destructively 
interferes. This is mitigated somewhat by the charged Higgs-top contribution, which always 
interferes constructively with the SM matrix element. The interplay between these two 
dominant SUSY contributions still leads to a predicted value below the observed one in 
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Figure 7: The statistical pull of some important constraints on mGMSB, where no dark 
matter constraint is applied. Posterior PDFs are plotted for: red (solid) lines have linear 
priors, green (dashed) lines have log priors, blue (dotted) lines have natural priors and cyan 
(dash-dotted) is the experimental likelihood constraint. The figures show (a) the top mass 
rrit, (b) the W mass rriw, (c) the logarithm of the branching ratio for Bg — * fi^fi^ with the 
black arrow being the current 95% CL experimental upper limit, (d) the anomalous magnetic 
moment (5a^, (e) BR{B — ^ Xg'y), (f) the branching ratio BR{Bu — * tu) divided by its SM 
prediction and (g) the isospin asymmetry in i? — ^ K*'y decays Aq-. 
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mGMSB. For natural priors the PDF is quite sharply peaked around the SM value of 
3.15 X 10~^. For log and linear priors the peak is shifted slightly to lower values, with the 
tails of the distributions being longer than that of the natural priors. 

Fig. [71(f) shows the posterior PDF of the ratio Rbtu = BR{Bu tv) / BR{Bu — ^ tu)sm, 
where BR{Bu tv)sm is the SM prediction of the branching ratio. The expressions for the 
SUSY corrections to this quantity can be found in [611 [62]. The SUSY contribution is always 
negative, so that the total MSSM prediction is smaller than that predicted by the Standard 
Model. The experimental value of R^g^^j = 1.28 ± 0.40 is compatible with a non-zero SUSY 
contribution. However, a more precise experimental measurement of this value which had 
the same central value but half the uncertainty would be better fit without SUSY, and could 
lead to some tension with the {g — 2)^ measurement. The isospin asymmetry in 5 — >■ K*"^ 
decays Aq- is shown in Fig. Ws)- mGMSB prefers a larger asymmetry than is observed. 
The discrepancy between the mode value and the experimental central value is just under 
2a. 

We have also investigated the PDFs for the electroweak observables calculated using 
SusyPope. We have found that they do not exhibit enough variation over the parameter 
space to present large constraints on our fits. In particular _R°, _R°, Ab, Ac and /Pj^ are 
all consistent with experiment independently of the prior. The forward backward asymmetry 
is predicted too large by 2.4 a while the left-right asymmetry is too small by 2 cr. 

To summarize the mGMSB parameter space fits, we observe a similar significant level of 
prior dependence in mGMSB to that previously observed in the CMSSM [3llll[71[6l[T6 l[T0l[TT] . 
It is clear from Figs. [3] and [H that none of the parameters are constrained independently of 
prior and so the fits require additional more precise and/or direct data before they can be 
considered robust. 

6.2.2 mAMSB 

Fig. [8] shows the ID posterior PDFs for some of the observables which constrain the mAMSB 
likelihood for the asymmetric Cdm- A comparison with Fig. [7| shows that many of the pos- 
teriors of the observables are similar in the mAMSB and mGMSB cases. Only and 
BR{B Xs'~f) show an large difference. The posterior PDF for rrit in Fig. [H](a) shows that, 
similarly to mGMSB, mAMSB is also in approximate agreement with the experimental con- 
straint. Fig. [8](b) shows the posterior PDF for mw which is under-predicted by mAMSB 
similar to the other models we have considered. Fig. M^c) shows the logarithm of the branch- 
ing ratio for the process B^ fi^fi~. The central SM prediction lies at the peak of each 
posterior. Our prediction neglects flavour mixing in the squark sector which result from the 
AMSB soft-breaking terms [26]. Ref. has shown that the effects of including squark 
flavour mixing can decrease BR{Bs —>■ /^"""/i^) by up to a factor of two for tan (3 > 22. This 
would result in a shift of logiQ{BR{Bs /U^/U") in the plot by -0.3. The 95% upper experi- 
mental limit placed on this process is < 5.8 x 10~® [65], which is always comfortably evaded 
even without this additional negative correction. 

The SUSY contribution to the anomalous magnetic moment of the muon is presented in 
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Figure 8: The statistical pull of some important constraints on mAMSB when the asymmetric 
dark matter constraint is applied. Posterior PDFs are plotted for: red (solid) lines have linear 
priors, green (dashed) lines have log priors, blue (dotted) lines have natural priors and cyan 
(dash-dotted) is the experimental likelihood constraint. The figures show (a) the top mass 
rrit, (b) the W mass rriw, (c) the logarithm of the branching ratio for Bg — * ^^^^ with the 
black arrow being the current 95% CL experimental upper limit, (d) the anomalous magnetic 
moment (5a^, (e) BR{B — ^ Xg'y), (f) the branching ratio BR{Bu — * tu) divided by its SM 
prediction and (g) the isospin asymmetry in i? — ^ K*'y decays Aq-. 



25 



Fig.[8]^d). With natural priors, the SUSY contribution to the anomalous magnetic moment is 
very small due to heavy sparticles and low value of tan j3. Light sparticles and intermediate 
values of tan/5 favoured by the log priors allow a better fit to experiment. Larger values of 
tan/5, while being consistent with 5a^, would then lead to a worse fit to BR{B — > Xs'j). 

Squark flavour mixing corrections, which we neglected, can in principle have an effect on 
the BR{B Xs'j) prediction in mAMSB, but the effect is small I26]. The BR{B Xs'y) 
posterior PDF shown in Fig. [8](e) approximately follows the experimental constraint for log 
priors, in contrast to mGMSB. The natural priors case follows the experimental constraint 
less closely and the linear priors somewhere in between the natural and log prior cases. This 
behaviour is consistent with the values of the evidence in Table IH where the log priors allow 
a better fit than the other two priors. 

The two observables shown in Fig. [8](f,g) are the ratio Rbtu and the isospin asymmetry 
in B meson decays Ao_, respectively. Both of these observables show the same behaviour 
as in mGMSB: not very constraining, due to the weakness of the large uncertainty in the 
SM prediction. While there is a small disagreement between the mAMSB prediction of Aq- 
and the experimental constraint there is too little variation of this observable over parameter 
space for it to mould the posterior significantly. The behaviour of the electroweak observables 
is qualitatively the same as in mGMSB, with all parameters except A^jf and A^j^ being in 
good agreement with the data while showing little variation over parameter space. As the 
accuracy with which these parameters are known will not improve in the near future, we are 
skeptical that they will prove useful in constraining the parameter space of simple models 
such as mAMSB or mGMSB. 

6.3 Combined 6a^—BR{B Xg^) constraint 

It is interesting that mAMSB is the only model we study for which the SUSY contributions 
to {g — 2)^ and BR{B X^'y) are both probably of the right sign and magnitude to satisfy 
the experimental constraints simultaneously. There are a number of factors which contribute 
to this. 

While the experimental and SM predictions for BR{B Xg'y) are consistent with each 
other, the fact that the experimental result is larger than that predicted by the SM means 
that given a SUSY contribution to this observable of a specific magnitude, a positive contri- 
bution will always agree with the data better than a negative one. mAMSB with yU > can 
accommodate such a contribution due to the sign of the coupling At. In mGMSB, mSUGRA 
and the LVS this will not be the case in the majority of parameter space. However, we note 
that using a more recent evaluation [86] of the Standard Model prediction which gives the 
branching ratio as 3.28 ± 0.25 x 10~^ would lessen this effect somewhat. 

A further effect is the hierarchies present in mAMSB spectra compared to the relatively 
compressed mGMSB spectra. Recall from Eq. [T3] that the magnitude of the supersymmet- 
ric one-loop g — 2 correction scales as ~ Mi^2^J'isL.ia(3/rn^, where m is the mass scale of the 
relevant particles in the loop, i.e. the chargino/sneutrino or neutralino/smuon mass scale. 
The dominant SUSY contributions to BR{B —>■ Xs'j) come from top-charged Higgs and 
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Figure 9: The 2D joint posterior in the {g — 2)^-BR{B X^'y) plane for (a) niAMSB and 
(b) mGMSB. Log priors and no dark matter constraint apply in both cases. 68% and 95% 
Bayesian credibility regions are shown as the inner and outer curves, respectively. The ellipse 
is centered on the experimental values of the two observables and show the 68% confidence 
level bounds on 6a^ and BR{B — > Xg'j). 



stop-chargino contributions. Agreement with the data typically prefers that the SUSY con- 
tribution to the branching ratio is not too large, i.e. the charged Higgs and stop masses are 
not too small. Thus, models in which the ratio of these masses to the slepton/neutralino 
masses is large have the potential to fit both {g — 2)^ and BR{B — > Xg'-f) simultaneously. In 
mGMSB, the spectrum is rather compressed compared to mAMSB, and so models with light 
neutralinos and sleptons also tend to have relatively light charged Higgs and stops. Thus, 
mGMSB fits the combination of the two observables less well than mAMSB. The interplay 
between {g — 2)^ and BRiB Xs7) fits has recently been explored in detail in [TT] for the 
case of the CMSSM, and previously also in [871 188] . 

We illustrate this behaviour in Fig. [9] which shows the 2D posterior PDFs in the ((7 — 
2)fj^-BR{B Xs'y) plane for mAMSB and mGMSB with log priors and no dark matter 
constraint. It shows the correlation between the two observables for mAMSB and anti- 
correlation for mGMSB, along with red ellipses showing the combined experimental con- 
straint. This (anti-)correlation is expected from the signs of At and M3 in the models |26j . 
Parts of mAMSB parameter space that provide a good fit to all of the data are also within 
the error ellipse, contrary to the mGMSB case. 

6.4 Best-fit points 

Fig. [ini shows the best-fit points (defined to be the highest likelihood point) for mAMSB, 
mGMSB, the CMSSM and LVS. The asymmetric DM constraint was taken, except for 
mGMSB where no DM constraint was applied. The caption contains the parameters of 
each model that yielded the best-fit, and the bars on the right-hand side measure how far 
each observable is from its experimental central value. We should note that our fits contain 
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Figure 10: Best fit points. The best-fit point prediction and experimentally measured values 
for each observable are listed. The pulls, represented by the horizontal bars, are in units of 
standard deviations. The corresponding parameter points are: mo = 312.7 GeV, 771^/2 = 45 
TeV, tan/5 = 15.9 for mAMSB, N^ess = 7, M^ess = 1-23 x 10^ GeV,A = 19.5TeV tan/5 = 
15.9 for mGMSB, mi/2 = 328 GeV, tan/5 = 11.6 for LVS and mo = 3338 GeV, mi/2 = 382 
GeV, Aq = 635, tan /5 = 8.6 for the CMSSM. fi is positive for all of the best-fit points. 
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quasi-degenerate likelihood maxima, and this combined with the fact that the MultiNest 
algorithm is not optimized to find the best-fit point, may mean that the parameters of the 
best-fit point are poorly determined. This is not important for any of the Bayesian inferences 
we have presented, but it is instructive to examine the properties of some point in each model. 
The point sampled with the highest likelihood is a good candidate for the point to choose. 
The CMSSM best-fit point is in the focus point region [Hj. It is not a really good fit to the 
anomalous magnetic moment of the muon, since heavy sleptons and moderate tan (3 render 
the SUSY contribution small. The other models' best-fit points all fit (5a^ much better, 
owing to their lighter spectra. Comparing the pulls of the electroweak observables, we see 
that there is not much variation in the values coming from them across the models except 
for M]y. The forward-backward asymmetry of e~^e^ bb, , shows a small change, being 
a little larger in the LVS than the other models. The pull of the left-right asymmetry A^p^ 
also shows some small dependence. All of the models listed have a predicted relic density 
of dark matter much smaller than the WMAP-inferred central value, requiring a dominant 
component of non neutralino dark matter (except for mGMSB, where we do not apply the 
DM constraint). 



7 Conclusion 

We ask the question: are any of the most commonly assumed low-parameter SUSY breaking 
mediation mechanisms favored over the others by current indirect and cosmological data? 
Aside from direct searches, SUSY corrections contribute to the observables we consider only 
in loop effects. For observables that agree with the standard model prediction, arbitrarily 
heavy sparticles suppress such loop effects and so will fit that observable. However, the 
anomalous magnetic moment of the muon prefers a significant contribution from supersym- 
metric loops, and so this observable prefers light SUSY particles. The dark matter constraint 
from cosmological observations cannot be accounted for by the Standard Model, but can be 
fitted by the MSSM if the lightest supersymmetric particle is the neutralino. We argue that 
the Bayesian evidence is the relevant statistical quantity for such an analysis. 

When performing such a statistical global fit, it is important to check robustness. The 
fits should be dominated by the data and not by the form of the prior if we are to claim 
robustness. Previous studies showed that fits to the CMSSM with four extra parameters are 
prior dependent [U [6l [10] , whereas the LVS model [13] (with only two extra parameters) is 
more robust. Thus it was natural to select mGMSB and mAMSB, each of which have three 
extra parameters to check robustness of the fits and compare the models against each other. 
One expects the level of robustness to go down with higher numbers of parameters, which a 
recent fit to the phenomenological MSSM [14J with twenty extra parameters, illustratecQ. We 
have presented constraints on the mGMSB and mAMSB parameter spaces, and found that 
in both cases there is significant prior dependence. Parameter inference from the models are 
therefore not robust and therefore require further more precise and direct data, perhaps from 
collider measurements of SUSY particles. Unsurprisingly, no robust statement regarding the 

"^There were a couple of prior independent inferences in the fits, such as the hghtest CP-even Higgs mass. 
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sign of fi can be made for mAMSB or niGMSB. We therefore encourage that future work 
on mAMSB and mGMSB include both signs of /i until the data is strong enough to support 
a prior-independent pronouncement on the status of the sign of /i, and not to disregard 
the /i < branch of the theory based on the preference of one observable for n > 0. The 
Large Volume Scenario, with two fewer parameters than the CMSSM is more constrained 
and shows a robust moderate preference for yU > 0. 

The model preferred by the data depends on what we assume for the DM relic density: 
whether it is made entirely of neutralinos (symmetric constraint) or whether we allow for 
the presence of non-neutralino dark matter (asymmetric constraint). An analysis of the 
constraining power of the various observables showed that it resides dominantly in the DM 
constraint in the case of the CMSSM and the LVS. This is not the case in mAMSB where 
the relic density is uniformly too small by an order of magnitude across parameter space, 
and the main constraint comes from the combined electroweak observables. Dropping the 
DM constraint altogether allows a comparison with mGMSB, but then no strong robust 
preference for any model can be found. However, for the symmetric constraint, mAMSB 
is strongly disfavoured (since it predicts essentially no neutralino dark matter) over the 
CMSSM and LVS. With the asymmetric constraint, mAMSB is at least moderately favoured 
over the CMSSM. 

Experience and familiarity with the methods of model selection and Bayesian inference 
from work such as that contained here will be invaluable once further more constraining data 
become available, hopefully from SUSY signals at colliders. 
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A Dark Matter Direct Detection 

In this appendix we present results on DM direct detection cross-sections for the Large Vol- 
ume Scenario. Direct detection rates in the CMSSM have been most recently been presented 
in Bayesian global fits in Refs. |H1 Ej. Our updated CMSSM fits are similar to the results of 
those articles, so we do not include them here. Our calculations of the detection cross-section 
have been obtained with micrOMEGAS2 . 3 . 1 [i7lli8| HOllSO] . 
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Figure 11: Posterior PDFs with linear priors for (a) the spin-independent cross-section a si 
and (b) the spin-dependent cross-section asD for neutrahno-nucleon elastic scattering plotted 
against the neutralino mass in the Large Volume Scenario, marginalized over both signs 
of /i. Both plots have a symmetric dark matter constraint. 

Neutralinos are Majorana fermions and can therefore interact with quarks via scalar 
and axial-vector interactions, but not via a vector-like interaction. In the scalar (spin- 
independent) case, the neutralino can interact at tree-level with quarks in the nucleus 
via Higgs or squark exchange. At one-loop the neutralino can couple to gluons via a 
quark/squark loop. At zero momentum transfer, the total cross-section is given by 



where rrired is the nucleon-neutralino reduced mass, Z and {A — Z) are the number of protons 
and neutrons in the nucleus and is the coupling of the neutralino to the proton (neutron) 
respectively. 

The neutralino-quark axial-vector interaction leads to a spin-dependent coupling propor- 
tional to J(J + 1), where J is the spin of the nucleus. The cross-section is given by 



where v is the relative velocity of the neutralino and T(f ^) is the scattering matrix element. 
This expression is then integrated over the Boltzmann velocity distribution of the neutralinos 
to obtain an average cross section. 

Since the results of our fits indicate that the LVS is sufficiently constrained so as to 
be independent of the prior, we show only the case of linear priors. The cross-sections are 
independent of whether the nucleon is a proton or neutron. We therefore have plotted the 
average (cr^ + cr"')/2 for each point sampled. As elsewhere in this article we marginalize over 
the sign of yU. 

Fig. [TT](a) shows the 2D posterior PDF for the logarithm of the spin-independent cross- 
section per nucleon in picobarns against the neutralino mass m-^. We also show the CDMS ^89j 
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upper bound assuming the local DM density to be po = 0.3Gev/cm^. Fig. [TT](b) shows 
the posterior PDF for the log of the spin-dependent cross-section per nucleon versus the 
neutralino mass. The DM detection cross-sections are both well below the bounds set by 
the current generation of DM searches such as CDMS^89j in the spin- independent and Su- 
perK [on] in the spin-dependent cases. Future one-tonne detectors should be able to probe 
the entire range of the LVS parameter space for the spin-independent cross-section down to 
(7si ~ 10~^°pb. The spin-dependent result is independent of the sign of p giving a robust 
prediction of gsd ~ 10"^-^=^°-^pb. 
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